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Abstract 
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photon-photon anomaly vertex deduced from available experimental data. 
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I. INTRODUCTION 



New measurement of the anomalous magnetic moment of the muon (muon anomaly) 
(ifi = I (5'/^ ~ 2) is underway at Brookhaven National Laboratory The anticipated level 
of precision, 40 x 10^^^, is more than 20 times higher than that of the best previous result 
i 

a^(exp) = 1 165 923 (8.5) x 10"^ (1.1) 

where the numerals in parentheses represent the uncertainties in the last digits of the mea- 
sured value. The primary purpose of the new muon g — 2 experiment is to verify the presence 
of the electroweak contribution. Other effects of potential interest are those of supersym- 
metric particles and leptoquarks 

At present the standard model prediction of is 

a^(th) = 116 591 714 (96) x IQ-^K (1.2) 

This consists of five parts: 

(i) Pure QED contribution Q 

a^(QED) = 116 584 705.7 (1.9) x 10~". (1.3) 

(ii) Hadronic vacuum polarization contribution obtained mainly from the measured hadron 
production cross section in e+e~ collisions PHT^- We quote here the latest evaluation 



that includes additional information obtained from the analysis of hadronic tau decay data 

mm 



^ Eq. (|l]^) is obtained from the measured value of the electron anomaly Ue 0] minus a small 
correction to Og due to muon, hadron, and weak interactions ||8|, plus the terms of dependent on 
the electron and tau masses evaluated using the fine structure constant obtained from the electron 
anomaly |p. 
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a^(had.v.p.) = 7 Oil (94) x 10 
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(1.4) 



(iii) Higher order hadronic vacuum polarization effect [|T^: 



(higher had.v.p.) 



101 (6) X 10 
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(1.5) 



(iv) Hadronic hght-by-hght scattering contribution HlBj : 



a^(had. 1-1) 



52 (18) X 10 
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(1.6) 



and a similar result obtained independently in Ref . [JT^ . 

(v) Electroweak contribution of one ||T^ and two |[TP| loop orders: 



a^iweak) = 151 (4) x 10 
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(1.7) 



The current uncertainty in the theoretical value of comes mostly from the hadronic 
contribution. It must be improved by at least a factor of two before we can extract useful 
physical information from the new high precision measurement and impose strong constraints 
on various candidates for possible extension of the standard model. 

The hadronic contribution appears for the first time in the order as the effect of 
hadronic vacuum polarization. (See Fig. 1 of Ref. for the Feynman graphs responsible 
for such a contribution.) Fortunately, the contribution of this type does not require ex- 
plicit ab initio calculation based on QCD, since it is precisely calculable from the measured 
hadron production cross section in e^e~ collisions [p|-[TT||. Future measurements at VEPP- 
2M, VEPP-4M, DA$NE and BEPS, as well as analysis of the hadronic tau decay data are 
expected to reduce the uncertainty of this contribution to the level that satisfies our need 



The contribution of the hadronic light-by-light scattering subdiagram is much smaller 
but is potentially a source of more serious problem because it has been difficult to express 
it in terms of experimentally accessible observables. At present it depends entirely on the- 
oretical consideration. This contribution has been estimated recently by two groups within 
the framework of chiral perturbation theory and the 1/Nc expansion [|l^,|l^. The leading 



pa- 
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terms arise from three types of diagrams shown in Fig. |l]: (a) pion-loop contribution, (b) 
pseudoscalar pole contribution, and (c) quark-loop contribution. The results obtained in 
Ref. [|16| for these diagrams are 

«m(«)Ihks = -4-5 (8.1) X 10-", 
«/^WIhks =-57.5 (11.4) X 10-", 

«m(c)Ihks = 9.7 (11.1) X 10-". (1.8) 



They add up to ( |1.6|) . A small axial- vector contribution 



(axial- vector) I = —1.74 x 10 , (1.9) 



was also obtained but not included in ( |1.6|) . ( See eq. (4.32) of Ref. ||T6|.) The corresponding 



results obtained in |jT^ were 



«m(«)Ibpp = -19 (13) X 10-", 
a^(&)lBPP = -85 (13) X 10-", 

a^(c)|Bpp = 21 (3) X 10-". (1.10) 



The effects of axial- vector and scalar poles were also considered in Ref. |117| : 



(axial- vector) I gpp = —2.5 (1.0) x 10 ", 
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(scalar) I Bpp = -6.8 (2.0) x 10-". (1.11) 
Summing up (|1.10| ) and (|1.11|) , they obtained 



a^(had. 1-1)Ibpp = -92 (32) x 10-". (1.12) 

The reasons for the difference between ( |1.8|) and ( |1.10|) are as follows: For the pion-loop 
contribution a^{a) it is due to the fact that the effective Lagrangian responsible for ( |1.10| ) has 
a non-derivative p°p°7r+7r- coupling while the corresponding term is absent in (|1.8D, which 



is obtained from a vector meson dominance (VMD) model with hidden local symmetry pi . 
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The absence of this term in the latter is a direct consequence of the fact that it satisfies the 
Ward-Takahashi identity [|1^] and a soft pion theorem for V^tt scattering amphtude . 



In Ref. [0] an effective chiral Lagrangian was proposed which reproduces the earher 
VMD result of Ref. |p in terms of an interaction term consisting of an infinite series of 
derivatives of the pseudoscalar meson. These higher derivative terms are accompanied by a 
mass scale Mc, which, after their resummation, becomes just a pole mass corresponding to 
a vectorial degree of freedom. Their result will be justified if that Mc can be identified with 
the vector meson mass Mp. However, the analysis of the low energy behavior of the p^n^ 
scattering prevents us from interpreting the above degree of freedom as being associated 



with the physical vector meson ||2^. The scale Mc will have to be associated with another 
physical degree of freedom, which is heavier than Mp. For such an Mc the magnitude of the 
charged pion loop contribution based on their model will become smaller than that given 



m 



17]. In this model the relation between the scale Mc and the resulting value of the 



charged pion loop contribution a^(a)lBpp in ( |1.1CI| ) is not transparent. For this reason we 
will henceforth choose ap{a)\^^^ of ( |1.8| ) as the contribution of Fig. |I| (a). 

The difference in the evaluation of a^(6) is mainly due to the fact that only the contribu- 
tions of the 7r° pole and rj pole are taken into account in ( p..8D whereas ( |1.1UD includes also 
the contribution of 77', which turned out to be not negligible When the 77' contribution 
is added to (|1.8|) (see Sec. the remaining difference for a^(6) is no longer large and 

reflects mainly the ambiguity and difficulty in carrying out the chiral perturbation theory 
estimate beyond the momentum range of several hundred MeV. 

Similarly, the difference between a^(c) of ( p..8|) and ( |1.10|) originates from the difficulty 
in estimating the contribution from large momentum region. 

As is seen from ( |1.8| ) and ( p..lO| ), the most important contribution comes from the dia- 
grams of type (b) in which neutral pseudoscalar mesons P { P = 7r°, rj and rj' ) propagate 
between two P77 vertices, as shown in Fig. 0. As is well known, the chiral anomaly for ttq, 
for instance, can be expressed by the effective interaction 
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-TT 



fJ.uXcr p p 



;i.i3) 



where /tt — 93 MeV is the pion decay constant, in the lowest order of chiral expansion. 
When apphed to the calculation of a^{b), however, this Lagrangian leads to an ultraviolet- 
divergent result. This divergence arises from the triangular photon-photon- muon loop which 
is obtained by reducing the 7r°77 vertex, represented by a large shaded blob on the left-hand 
side of Fig. |^(a), to a point. (Fig. ^(b) is convergent in the same limit.) It is a signal 
that the local interaction ( 1.13| ) is not applicable to photons and pions far off mass shell. In 
fact, such a triangular diagram will have a damping behavior in the underlying QCD theory, 
which protects its contribution to the muon g — 2 from diverging. 



In Ref. 1 16], four models were considered to examine the effect of various assumptions 
on the off-shell behavior of the 71^77 form factor: 

(61) vector meson dominance (VMD) model, 

(62) quark triangular loop (QTL) model, 

(63) QTL model combined with the VMD model. 



(64) extended Nambu-Jona-Lasinio (ENJL) model p3 



The contribution a^(6;7r'') of Fig. ^evaluated by these models, and given by Eqs. (4.1), 



(4.3) and (4.4), and Table V (with gA = 0.5, Mc = 1 GeV) of Ref. |]T6l, are reproduced here: 



-55.60 X 10-^\ 
-86.90 X 10"", 
-33.76 X 10-", 
-42.84 X 10"". 



:i.l4) 



Similar calculations were carried out in Ref. for the t] resonance, too. The final value of 
a^{h) reported in Ref. |TB| and quoted in (Of) is based mainly on the model (64). 



^The VMD model can be justified within the hidden local symmetry picture of the chiral model 
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After completion of these calculations an experimental measurement of the P77* form 
factor, where P stands for 7r°, rj or r/', came to our attention [p^-p6|. This is very important 
in the sense that it opens up the possibility of evaluating a^{h) utilizing the experimental 
information instead of relying solely on theoretical consideration. This may enable us to 
reduce significantly the uncertainty in the evaluation of a^{b). 

The purpose of this paper is to amplify the preliminary discussion in Ref. about the 



implication of the measurements and discuss in full the effect of the measured P77* form 
factor (Sec. |T|) and possible impacts of the yet-to-be-measured P7*7* form factor on the 
evaluation of a^(6) (Sec. |rv|) . Before going into these sections, we review the theoretical 
aspects of the asymptotic behavior of the form factor of our interest in Sec. |T[ The last 
section (Sec. 0) will be devoted to the summary and discussion of our results. 

II. THEORETICAL OFF-SHELL tt^j*j* FORM FACTOR 

Let us write the invariant P7*7* amplitude as 

-M(7*(Pi,Ai)7*(P2,A2)-P(g)) 

= 4MyxM)e,uapp'ip^2Mp{pl pI q' = ml), (2.1) 

and define the form factor Fp{p\, pi, q^) by 

Fp{pl pI q') = ^Mp{pl, pI, q'). (2.2) 
a 

In the chiral limit (g^ = m?p = 0), this is normalized as 

Pp(0, 0, 0) = 1. (2.3) 

Let us first examine the off-shell behavior of the 7r''7*7* form factor in the quark triangular 
loop (QTL) model. In this model, we find (ignoring isospin violation) 

FS^^pI, pI, q') = I^iipl, pI, q') 

[d^]— -2 2^ (2-4) 

mi - Z2Z3pi - zszipi - ZiZ2q^ 



where [dz] = dzidz2dz^5{l — zi — Z2 — z^), and tjIu ~ 300 MeV is the constituent mass of up 
(and down) quark ||2^. For pi = pi = q"^ = this function reduces to (^.31). Carrying out 
the integration one can readily find that 



FTiPl Pi 0) 



m 



Pi -pI 



In 



In 



Ami -Pi + ^-Vi 



Ami- Pi - y-Pi 



Ami ~ P2 



-pI 



Ami ~ P2 



-pI, 



(2.5) 



For large p\ with p\ = Q and = this has the asymptotic behavior of the form 



Fj"(p?, 0, 0) 



For large p\ ~ p\ and = we find 



-Pi 



In 



-pY 



m: 



(2.6) 



FriPlplO)r.H\n 



-pX 



-Pi 



mt 



(2.7) 



Note that both ( |2.6| ) and ( p.7| ) have the same power behavior at large momentum transfer, 
differing only in the logarithmic factor. 
Let us define the function jF^o by 

1 



Ana 
1 



MA-Q\ 0, 0) 
FA-Q\ 0, 0) 



(2.8) 



Fig. 1^ shows the momentum dependence of Q'^J-'-n-o^Q^) for the cases (61), (62) and (63) 
together with the experimental data reported by CLEO collaboration . This shows that 
the VMD model (61) fits the experimental data particularly well. The model (62) gives 
slower damping for higher momenta. This is because of the extra (ln(— p^))^ factor found in 
( p.6[ ) which compensates for the leading damping l/pl to some extent. The model (63), on 
the other hand, predicts rapid decrease due to the stronger damping factor (ln(— j9^))^/p^. 

Nonperturbative analysis of asymptotic behavior of the exact form factor F^o(p^, p^, q^) 
has also been carried out by several methods P^-pS. Various corrections to the result of 



have also been considered: Ref. [30| discusses the effect of the parton transverse momentum. 



Ref. [0 considers the gluonic radiative corrections to [^]. Ref. p2| derived an asymptotic 



formula applying the QCD sum rule |33|. These corrections have been compared with the 



data in Here, we concentrate on the methods of Refs. [p^^-|29[] since they are more 

directly related to the consideration of this paper. 

For P2 = q"^ = ^ all three methods [P7|-p5| predict the same leading power of momentum 
for large pf. But the coefficients found are different: 

3 2n'p^ 



hm FApi,0,0) 



-pI ' 
1 2vrV.^ 



[28] 



29 



(2.9) 
(2.10) 
(2.11) 



Nc -pI ■ 

The asymptotic behavior of F^o(p^, p^; 0) fo^ large p\ = pi has also been studied nonper 
turbatively ||34]| : 



1 27r2/,' 



lim F^o{pI,p10) 



(2.12) 



-Pl^oo Nc -pi 

Ref. states in addition that ( p.ll| ) holds for all values of p^ including p^ ^ p\. 

Different asymptotic behaviors of ( |2.9| ), ( |2.1CI| ) and ( |2.11| ) presumably reflect different 
physical assumptions. Eq. ( |2.9| ) is obtained by appealing to the parton picture in the 
infinite momentum frame. The leading momentum power dependence there comes from the 
scaling behavior of the pion wave function. However the coefficient depends on the long 
distance aspect of the pion, which is vulnerable to theoretical prejudices. According to 



literatures P6| , |30| , |31| the ansatz adopted by Brodsky and Lepage for the pion wave function 
is not badly in conflict with the experiment. 

Eq. ( p.lO| ) was obtained using the operator product expansion (OPE) technique, within 
the approximation of the dominance of operators of lowest dimension and twist The 
reliability of OPE and truncation of operators depend on the detailed prescription of how 
the large momentum limit is taken. In its application to the current issue the form factor is 
expanded in = — (pi +P2)^/4 and uj = {p\~p\) /Q"^ ■ The leading power and its coefficient 
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can be determined including long distance effect only when the absolute magnitude of u is 
small. This is the case for large p\ ~ p\ but not for large p\ with p\ fixed to 0. In the latter 
case the estimate based on the OPE is not trustworthy since uj will not be restricted to small 
values and the long distance contributions associated with higher power of uj might become 
as significant as the leading term P? . 



The derivation of (|2.11| ) relies on the Bjorken- Johnson- Low theorem This theo- 



rem is usually used to define the commutator of operators by postulating an appropriate 
asymptotic behavior for some correlation function. The asymptotic behavior in on 



the other hand, is derived from the Bjorken- Johnson-Low (BJL) limit of the matrix element 



of the commutator of the electromagnetic currents ||29|. The use of equal-time commutation 
relation of the quark fields turns this commutator into the axial current. This enabled them 
to express the relevant matrix element in terms of the pion decay constant f^^. However, 
this approach may not be reliable because the canonical equal-time commutation relation 
does not necessarily lead to the correct anomaly relation |^. Whether this caution is also 
relevant for the evaluation of the commutator of the electromagnetic currents or not will 
not be pursued here. Irrespective of the validity of derivation of ( |2.11| ) in [^, we note 
that the asymptotic behavior of 71^77* in the time-like region will differ significantly from 
that measured in the space-like region. At present we cannot find reason to support the 
possibility that such a non-trivial continuation occurs between the two regions. 

In order to compare the theoretical predictions of the 71^77* form factor with experiment, 
it is necessary to interpolate them from the asymptotic region to pi = 0, where they are 
fixed by the anomaly condition. Ref. proposed a one-parameter interpolation formula 



where = —pi- The formula ( p.l3| ) has the same dependence as if we choose 
Stt^/^ = M^. Thus this interpolation reproduces the dependence very well. However, 
this also fixes the value of jF^o(O) which differs from the observed value by a factor of ~ 1.2. 
In this sense (|2.13| ) does not fully reproduce the experimental data. Interpolations of (|2.10|) 
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and ( |2.11| ) will be even more difficult in this respect. 



III. P77* FORM FACTOR 



A. Experimental results 

Depicted in Fig. ^ are the function Q'^J^tto{Q'^), corresponding to the asymptotic form 
factors (|2.9|) , (|2.1CI|) and and the interpolation function ( p.l3|) , together with the 

observed 7r°77* form factor provided by the recent CLEO data [^. It is seen that none 
of the asymptotic behaviors shown there is fully consistent with the new data provided by 
CLEO for = 93 MeV. 



On the other hand, it was noted ||2^ that the experimental data can be represented very 
well by the empirical formula 

647r r(7r'^ 77) 1 



J^ttO {Q 



(3.1) 



\ {Anaymlo 1 + (QVA^o) ' 
if one chooses A^o ~ 0.77 GeV and r(7r° — * 77) = 7.78 (56) eV. The uncertainty in T is 
about ±7% Note that this A^o is nearly identical with the physical p mass. 

Similarly, the experimental data for the rj and 77' form factors can be represented well by 

647rr(P^77) 1 



(3.2) 



for P = T] and rj' , with the choice of A^ ~ 0.77 GeV and A^/ ^ 0.85 GeV. Here, to take 
account of non-zero mass of meson properly, we have modified ( p.8| ) as 

1 



Ana 



Mp{-Q\ 0, m'p) 



(3.3) 



The life-times chosen are T{ri — > 77) = 0.46 (4) keV, and r(?7' — > 77) = 4.26 (19) keV. [p7 



B. Effect of the P77* form factor on muon g — 2 



Let us now evaluate the pseudoscalar pole contribution to the muon g — 2, assuming 
that all virtual photons are modified by the P7*7* form factor in accord with the VMD 
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model. We consider two models. The first model (VMB^) assumes that the vector meson 
dominance is realized in terms of the p and mesons of the form 



with 



fp TT 



(3.4) 



(3.5) 



for V = p and 0, while i^,j){P) is treated as an adjustable parameter. Here we have adopted 



the approximation Mp ~ M^. Note that fp appearing in (|3^ ) is not the usual decay 
constant. Rather it should be regarded as an "effective" decay constant. When we fit ( p.4|) 
with the experimental data we use the expression 



^.o(Q2) 
which follows from (|3.3|) and 



647r r(P ^ 77) \Mp{-Q\ 0, ml) 

~ X 



^ (47ra;)2m 



|Mp(0, 0, ml) 



(3.6) 



r(p ^ 77) 



mp 



Mp(0, 0, ml) 



(3.7) 



Note that fitting of ( p.4|) with data does not involve fp. The value of fp itself will be 
determined once we fix k^{P) from the experiment for a given r(P — 77) through the 
relation 



r(p ^ 77) 



647r 



Mp{0,0,ml) 



1 mp 



647r Vvr 



a 



p;--(0, 0, m^) + «:,(P)Pr^-(0, 0, m^) 



(3.8) 



The second model (VMDp) assumes that the vector meson dominance is realized in 
terms of the experimental pole mass Ap obtained in |2^. We vary Ap within its error bars. 

In the model VMD^ (or VMD^) the parameter ^^(P) (or Ap) is determined by fitting to 
the experimental data for each pseudoscalar meson P and a fixed value of r(P ^ 77) 
using the CERN library routine MINUIT. Here the confidence level of la is imposed to 
extract the uncertainty associated with the parameter ^^(P) or Ap. We then examine the 
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effect of pseudoscalar poles on the muon anomaly based on this confidence level and the 



values of r(P — > 77) given in |]3^ 



The results of our calculation are summarized in Tables |— |VI[ They have all been 
obtained by the integration routine VEGAS with 800 thousands sampling points per iteration 
and iterated 15 times. From the explicit calculation, the dependence on r(P — 77) and 
K(^(P) (or Ap) is found to be rather simple, as can be partly inferred from those tables. To 
obtain the most likely value and the error of a^{P) we follow the procedure described below 
for the VMDa model. First we take the average of a^(P) for (r(P 77)max, «^fli(-P)ccnter) 
and that for (r(P — > 77)mm, ^^tt>{,P) center) and deduce the errors associated with r(P — > 
77) as the maximally allowed deviation from their mean value. Likewise the uncertainty 
associated with ^^(P) can be found from a^{P) for (r(P 77) center, i^<f>{P)msLx) and that 
for (r(P 77)center, i^<f>{P)mm) ■ The Central value of a^{P) is set equal to the mean value of 
the above two averages. The total uncertainty is obtained by taking the square root of the 
sum of the squares of the two errors, assuming that the errors associated with r(P ^ 77) 
and ^^(P) are independently distributed. The same procedure also applies to the VMD^ 
model with (r(P ^ 77), A(P)) instead of (r(P 77), k^(P)) . 

From the results listed in Tables §— we find that the contributions to the muon anomaly 
from 7r°, 77° and rj' propagation in the VMDyi model are given by 



a^(7r", VMDa) = -0.045 8 (29) x 



7r 



a^{ri, VMDa) = -0.010 6 (9) x (^^^ 



3 



a^{r]', VMDa) = -0.009 4 (7) x ) . (3.9) 

For the VMD^ model we find similarly 

a^(7r°, VMDb) = -0.045 8 (29) x ( 3 
a^{r], VMDb) = -0.010 7 (8) x 

a^(V, VMDb) = -0.009 5 (7) x (^^^^ (3.10) 
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The total contribution to muon g — 2 from the pseudoscalar pole effect becomes 

a^(6, VMDa) = -0.065 8 (32) x f 



TT / 

-11 



= -82.5 (4.1) X 10-'\ (3.11) 
and 



a^{b, VMDb) = -0.066 (31) x 



vr 



= -82.7 (3.9) X 10-^\ (3.12) 

where the errors above were deduced assuming that the uncertainties associated with dis- 
tinct mesons are normally distributed. These results are rather insensitive to the difference 
between the two models, both of which are already strongly constrained by the experiment. 

Let us comment here on the pion decay constant. Note that every postulated form factor 
used in the previous analysis is proportional to l//p implementing the PCAC Thus the con- 
tribution to from each pseudoscalar meson propagation becomes proportional to 1/fp- In 
the previous analysis P,p!B|, = 93 MeV was used. This is the decay constant of the charged 
pion. That of the neutral pion is quoted as fj^o ~ 84 (3) MeV in |^ . Correspondingly the 



values of predicted by all the above form factors receive the enhancement factor 1.22 if 
the latter value is used. But this implies that too much isospin breaking (about 10 %) shows 
up in the pion decay constants resulting from the mass difference of up and down quarks 
and electromagnetic corrections. Those values found in Refs. [0,^ for 7r° were extracted 



from the use of pole-type form factor ( |3.1| ) corresponding to the VMD picture together with 
A^o = Sn'^f^o. This relation must not be trusted once we treat T{P — > 77) and Ap as two 
independent parameters. Its application might lead to a fictitiously large isospin violation in 
the pion decay constant. In our analysis throughout this paper, the decay constants are not 
the direct input, because its determination primarily depend on the form factor assumed. 
Instead, we use the partial decay width r(P — 77) given in and determine the effective 
decay constant according to an equation such as (|3.8|). This provides a phenomenological 
and totally self-contained algorithm. 
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IV. EFFECT OF THE P7*7* FORM FACTOR ON MUON g-2 



For a full evaluation of the muon anomaly a^{h), we need information on the form factor 
for all values of pi and p2- In other words, we need information not only for large p\ 
with p\ = Q oi large p\ with p\ = Q but also for the case where both p\ and p\ are large 
simultaneously. In the VMD model (61) the off-shell 71^7*7* form factor for large p\ and p| 
takes the form 

^.o(Pi,P2,0) oc ^ for pIpI-p". (4.1) 

On the other hand the QTL model (62) or a nonperturbative estimate shows much slower 
behavior (~ 1/p^) as is seen from ( p.Tj) and (|2.12|) . 

This suggests that the 71^7*7* form factor may have hard components for large p\ and p\. 
Taking the asymptotic form ( p.l2| ) into account, we may, for instance, choose the function 

as a candidate of the hard component. Another possibility is to use ( |2.7| ) itself. With this 
in mind we consider two models: 



(I) Linear combination of the VMD^ model and the function 

(II) Linear combination of the VMD^ model and the QTL model (62)- 

We will not discuss the VMDb model since the difference between the VMDb model and 
VMD A model is small. Let us introduce a parameter (3p io describe the deviation from the 
VMD A model. The best value oi (3p depends on the species P of pseudoscalar meson. To 
be more specific, we write the scalar part of the invariant amplitude (|2.1| ) as 



Mp{pl pI q') = ^- [F^^pI pI q') + «:,(P)Ff (p?, pi q')] , (4.3) 

where, for the extended model (I), 

F^^ipl pI q') = (1 -/5P)F™(P?, P^) +/5PFS(P)(P?, pI), (4.4) 
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and, for the extended model (II), 

F^^ipl, pI q') = (1 - /5p)Fr°(p?, pI) + f3p F^^\pI pI q'). (4.5) 



Here F^^^{pl, pi) is the form factor (U) in VMB^ model, F];}l^{pl, pfj is given by ( ^ 
in which Mx{P) is assumed, for simplicity, to take the same value for p and 0, and 

Fv [Pi, P2, Q ) = { (4-6) 

Imsipl, PI, q^) i0TV = (j). 



Irriu is defined in ( |2.4| ). takes account of the fact that meson is dominated by the ss 
state. Again the "effective" decay constant fp in ( [4.3|) does not have the usual meaning and 
is determined in the same way as in Sec. PIB| . In the extended model (II), for instance, 
for each meson P and a given value of r(P —>■ 77), the direct fit to the experimental data 
determines the parameters Pp and K(f,{P), which give the least value of according to ( p. 61) . 
Once the values of T{P — > 77), Hct,{P) and Pp are fixed, a relation similar to ( p.8|) leads to 
fp- 

Note that the extended model (I) of ( ^.4| ) has an incorrect asymptotic behavior for finite 
Pi or p^. Thus it must be regarded as an empirical formula designed to fit the experimental 
data and applicable only up to the largest momentum transfer of the experimental data. 
We use it nevertheless because of simplicity and because its effect on the muon g — 2 is 
small as far as /5p is small. The extended model (II) is also not very satisfactory since the 
QTL part has a logarithmic growth which distorts the form factor even for moderately large 
momentum transfers. Thus both must be treated with some caution. 



A. Extended model (I) 

The extended model (I) involves one additional parameter Mx, and the best fitting values 
of i^^{P) and (3p change as Mx changes even for the same value of V{P — > 77). To avoid 
time-consuming and excessive analysis, we first perform two parameter fitting for fixed Mx 
and r(P ^ 77) to determine a set (/5p, k^{P)) that gives the best x^. Next we determine 
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the range of allowed at Icr confidence level with f3p fixed to the best value obtained above, 
and vice versa. The result of such a fitting is summarized in Tables |V11| — PC[ They show 
that (3 pis small in general and thus does not favor strong presence of hard component. Note 
that a large Mx implies near g-independence of the term: 

^ 1 for Mi ^ oo with fixed q^. (4.7) 



Mi -g2 

This means that the formula ([4.4|) tends to a constant (3p and hence gives rise to a divergence 
of the muon g — 2 proportional to InMjf. Larger Mx also leads to larger minimum of x^. 
Thus we restrict Mx to the values below 6 GeV, and examine the effect of this term on the 
muon g — 2. 

The computation of muon g — 2 has been performed for several sets of (/?, k,^) , in which 
either P or is its best fitted value for given Mx and T{P — > 77). The results listed 
in Tables |X|— pCII| for the extended model (I) have been obtained by integration with 800 
thousands sampling points per iteration which is iterated 15 times. The errors generated by 
the numerical integration itself are not explicitly stated since they are all under 0.05 %, far 
below the errors arising from the uncertainty in r(P — » 77). 

Let us first focus our attention on the case Mx = Mp in Tables |X|— pCll| , which is the 



lowest value of the scale Mx characterizing the modification of the high energy behavior in 
the present model. To deal with these data we adopt the same procedure as in the case of the 
VMD model. For instance the uncertainty associated with the variation of r(P 77) can be 
read off from a^(P) for (r(P 77)max, /5(^)center, n<t,iP) center) and that for (r(P 77)mm, 
/3(-P)center, «^fli(-P)ccntcr) • The Uncertainty associated with the variation of /3(P) and ^^(P) 
will be found similarly. We take the average of three tentative central values to obtain the 
central value of a^{P) given below. The combined error is obtained by taking the square 
root of the sum of squares of the errors. The contributions to muon anomaly from 7r°, 77° 
and rj' propagation in the extended model (I) can then be written as 



a^(7r°, (I)) = -0.045 6 (28) x 
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a,ir], (I)) = -0.010 3 (9) x (^^) , 

a^(r^', (I)) = -0.009 (6) X (4.8) 

where the uncertainties are estimated assuming that the error distribution associated with 
r(P — > 77), Pp and Kfj,{P) are normal and independent. The result obtained assuming Mx 
= 2 GeV is nearly identical with ( [4.8| ). The result for larger values of Mx deviates only 
slightly from ( |4.8| ), except for the considerable Mx dependence seen in a^{ri, (I)). 



B. Extended model (II) 

The trial form factor for the extended model (II) is the VMD^ model augmented slightly 
by a QTL term. The additional term receives a strong constraint from the experimental 
data as is shown in Tables PCIII| — pCV[ This is anticipated because the logarithmic factor in 



( p.7|) of the QTL model is quite visible even in the moderately large momentum region, and 
hence receives a strong constraint from experiment. Tables PCV1| — pTVlIl] show the prediction 
for the muon g — 2 from the trial form factor ( |4.5| ) for each pseudoscalar meson under the 
constraint of the parameters in Tables pCIII| — [XV| , in which all the results have been obtained 
by the integration with 500 thousands sampling points per iteration which is iterated 25 
times. In the manner similar to that of Sec. [IV A| , the contributions to muon anomaly from 
7r°, rj^ and 77' propagation in the extended model (II) are estimated to be 



a^{n", (II)) = -0.045 6 (28) x - 



a 



IT 



3 



a^{r], (II)) = -0.006 4 (9) x 



TT 



TT 



a^(V, (II)) = -0.008 7 (7) X - . (4.9) 



It is seen that the QTL modification to the VMD results of Sec. |ITT| is larger compared 
with the case of extended model (I). In particular, the t] contribution is reduced in magnitude 
compared to that obtained in the other scheme. Comparison of ( f4.9| ) with ( [4. 8] ) shows that 
the additional term in case (II) leads to a better x^- Moreover the signs of for rj and 



rl tend to become negative in the present case compared to the VMD^ model, as can be 
seen from Tables |I|, |T|, pClV| and pCV| . The opposite signs of p and contributions result 
in the decrease of absolute magnitudes of ri and v( contributions. This relative minus sign 
is required to fit to the experimental data once the form factor contains a QTL component 
which enhances the contribution of large momentum region. In the sense that this region 
is beyond control of the experimental data and the logarithmic enhancement may be an 
artifact of perturbation theory, however, the result ([4.9|) must be treated with some caution, 
and should be accorded less weight than the extended model (I). 

Note that the fiavor SU(3) relation implies /«0(?7) ~ —0.28 and ^^(r/') ^ 0.29 for d ~ 
—10.1° While the signs of K^^ij]) and ^^(V) iii the extended model (I) are consistent 
with the SU(3) relation, those in the extended model (II) are not. Of course this does not 
mean that one model is more appropriate than the other. In the first place, it is obscure 
whether SU(3) breaking effect due to strange quark mass is substantial or not. Secondly, 
we do not know the relevance of the U(l)^ anomaly (QCD anomaly) contribution to the 
r/Q — i> 77 (r/Q is the SU(3) singlet component) vertex. These questions have been examined 
for P — i> 77 in and for another decay channel in in which the subtlety that the 
presence of QCD anomaly raises is explored in identifying the decay constant invariant 
under renormalization group, inducing a nonperturbative ambiguity associated with each 
decay channel. We shall not discuss its relevance quantitatively here. We believe that these 
effects can be incorporated in a few constants parametrizing the form factor we chose, by 
fitting to the observed pseudoscalar production cross section through the two photon process 
2l-p6 . 



V. SUMMARY AND DISCUSSION 

In this paper we reexamined the pseudoscalar meson pole contributions to the muon g — 2 
taking account of the measured P77* form factor and of possible effect of the postulated 
P7*7* form factor. Insofar as we demand that these form factors respect the measured 



19 



strength of P77* for large momentum (~ several GeV) with one photon on the mass-shell, 
the form factor cannot deviate from the VMD model substantially. 

The results ( [4.8| ) and ( |4.9| ) of the extended models have somewhat smaller overall uncer- 



tainties than the results ( |3.11| ) and ( |3.12| ) of the VMDa and VMDb models. This may be 



due to better flexibility of the extended models. 

Inspection of ( [4.8|) and ( [4.9|) shows that the contribution of the 7r°7*7* form factor is 
nearly identical for both models while the contribution of the 777*7* form factor shows 
considerable dependence on the assumed structure of hard component. This may not be 
surprising since the effect of hard component will be stronger for rj than for n^. Better 
agreement of the ?7'7*7* cases may just be a matter of coincidence. The t]' contribution in 



( [4.8|) is somewhat smaller than that of rj, in contrast to the results of the ENJL model and 
the extended model (II), both of which contain the QTL form factor at least partly. 

Note that either ( |4.^ ) or ( |4.9| ) cannot be taken as our best estimate based on the exper- 
imental data only since they depend on theoretical assumption about the hard component. 
Instead they must be regarded as providing a measure of uncertainty in due to the un- 
known effect of the P'~f*'~f* form factor. Our calculation shows that this uncertainty comes 
mainly from the contribution of ri'y*'y* form factor. For the reasons discussed in Sec. [IV B| , 



we believe that a;x(?7, (II)) is less reliable than a^(r7, (I)). But their difference 0.0039 (a/vr)^ 
may be regarded as a measure of theoretical uncertainty. 

Based on these considerations, we adopt ( |3.12|) as our best estimate and choose as its 



uncertainty the statistical combination of the error in ( p.l2|) , which comes from the ex- 



perimental uncertainty in the P77* measurement, and the effect described above due to 
possible presence of hard component in the P7*7* form factor. We thus arrive at the total 
pseudoscalar pole contribution 

a^{n + ri + ri') = -82.7 (6.4) x 10"^^ (5.1) 

The uncertainty here has been reduced to one half of the previous result ( |1.8D . That the 
result ( |1.10|) from Ref. is close to our result (|5.1| ) is not surprising since the modified 
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ENJL model of [T^ contains an adjustable free parameter. The crucial difference between 
the old results ( |1.8| ), ( [1.10| ) and the new result ( |5.1| ) is that the latter is much less dependent 



on the theoretical ambiguity, being determined by the experimental data to a large extent. 

The theoretical uncertainty in ( |5.1| ) reflects the absence of data on the P'j*'-/* form factor. 
It is useful to note, however, that it is possible to gain some insight for the effect of the 
P7*7* form factor based on the measurements of 7r° e^e~ and r] — > decays, both 



of which involve the P'j*'j* form factor. The measured branching ratios ^ 



5,ncas(r/ ^ fi^l^-) = ^ ) = 1.4 (0.2) X 10-^ (5.2) 

and 

B^eM' ^ e+e-) = ^^^°o"tf^^)^ = '^•3 (1-9) >< 10~'' (5-3) 

where the latter is a weighted average of the data from Ref. [0] are in good agreement with 
the theoretical values based on the VMD model ^2| |^ 

5vMD(vr° ^ e+e') = ^."^''^''.^ = 6.41 (0.19) x 10~l (5.4) 

The differences between these values are presumably due to the effect of the P7*7* form 
factor. It appears to be no greater than 15 %. In principle, such an information can lead 
to a completely model-independent evaluation of the pseudoscalar pole contribution to the 
muon g — 2. Further improvement in these measurements will thus be of great interest. 

Unlike the pseudoscalar pole contribution discussed in this paper, the other contributions 
of hadronic light-by-light scattering type from Fig. ^ remain vulnerable to various theoret- 
ical ambiguities. The axial-vector pole contribution (|1 . 1 1|) was obtained entirely based on 



the ENJL model, which is certainly not satisfactory in the sense that it does not lead to 



the correct asymptotic behavior for the tt 77* vertex [|T7|. Based on the estimates (p..9| 



^We thank J. F. Donoghue for calling our attention to this reference. 
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and (|1 . 1 1|) the axial-vector pole contribution seems to be relatively minor compared with 
the effect ( |5.1| ) of the pseudoscalar pole propagation. Since nothing better is available at 
present and since it is very small, we choose ( p,.9| ) as our best estimate of the axial vector 
contribution assigning an uncertainty as large as the value ( |1.9| ) itself. The axial-vector 
meson contribution may be estimated more reliably if the off-shell structure of ^77 vertex 
is available from the experiment as was the case for the P77 vertex. 

Another possible pole effect, the scalar resonance inferred within the ENJL model, has 
not been observed convincingly. It may not be surprising if an "exact" QCD calculation 
leads only to a continuum spectrum instead of producing a broad resonance in the scalar 
channel of the corresponding energy scale ( 0.6 ^ ^/s ^ 1.5 GeV ). 

In the absence of any resonance in the low energy region (below several hundred MeV) in 
the TT — TT and qq channels, the pion-loop contribution a^(a) and the quark-loop contribution 
a^(c) may be regarded as representing the effect of the continuum spectrum. Previous 
analysis p6|jr7[] shows that a^{a) and a^(c) have sizable contributions from the region of 
higher momentum transfer where neither chiral perturbation theory nor the 1/Nc expansion 
may provide reliable guidance. Even mutual independence of the three types of diagrams of 
Fig. |I| might not hold valid there. For lack of better argument, however, we choose the sum 
of a^(a) and crude estimate of the entire continuum contribution. The scalar pole 

contribution may be regarded as included in this estimate. 

Finally, collecting a^(a) and a^(c) from (pT^ ), the pseudoscalar meson pole contribution 
from (|5.1|) , and including ( |1.9|) , we present 



a^(had 1-1) = -79.2 (15.4) x 10"'' (5.5) 

as our best estimate of the total hadronic light-by-light scattering contribution to the muon 
anomaly. Here the uncertainty is obtained from the errors of various components assuming 
normal distribution. If the uncertainties are combined additively, one would obtain ~ ±30 x 
10-11. 

Further progress in theory would have to wait for the lattice QCD calculation of the 
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hadronic four-point function. In view of the recent progress in the lattice QCD and rapidly 
increasing computing power, such a calculation may no longer be a far-fetched dream. As 
of now, the theoretical value of the muon g — 2, which consists of (|1.3| ), ( |1.4| ), (|1.5|), (|1.7D, 
and our new analysis of the hadronic light-by-light scattering effect ( |5.5| ), is given by 

a^(th) = 116 591 687 (96) x 10"^^ (5.6) 

The largest source of theoretical uncertainty remains to be the hadronic vacuuum polariza- 
tion effect ( p.. 41) . Better measurements of this effect are urgently needed. 
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TABLES 



TABLE I. 7r° contribution to muon 5 — 2 in the VMD^ model. 









7.22 


-0.01 


-0 041 87 f2) 


7.22 


0.12 


-0 043 56 f2~) 


7.22 


0.27 


-0.045 28 (2) 


7.78 


-0.07 


-0.044 03 (2) 


7.78 


0.04 


-0.045 80 (2) 


7.78 


0.16 


-0.047 54 (2) 


8.34 


-0.12 


-0.046 17 (2) 


8.34 


-0.03 


-0.047 95 (2) 


8.34 


0.08 


-0.049 76 (2) 


TABLE IL 


r/ contribution to muon — 2 in the VMD^ model. 






\ 7? ) 
"'<p V// 




0.42 


0.04 


—0 009 646 f4') 


0.42 


0.19 


-0 010 235 f4) 


0.42 


0.38 


-0 010 843 i4>i 


0.46 


-0.04 


-0.009 910 (4) 


0.46 


0.07 


-0.010 726 (4) 


0.46 


0.22 


-0.011 349 (4) 


0.50 


-0.10 


-0.010 59-1 (1) 


0.50 


-0.01 


-0.011 211 (4) 


0.50 


0.13 


-0.011 921 (5) 
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TABLE III. ri contribution to muon 5 — 2 in the VMD^ model. 



r(ry' ^ ^-y) [keV] 






4.07 


0.25 


-0.008 584 (4) 


4.07 


0.48 


-0.009 190 (4) 


4.07 


0.82 


-0.009 818 (5) 


4.26 


0.19 


-0.008 766 (4) 


4.26 


0.39 


-0.009 380 (4) 


4.26 


0.67 


-0.010 013 (5) 


4.45 


0.13 


-0.008 947 (4) 


4.45 


0.31 


-0.009 562 (4) 


4.45 


0.55 


-0.010 205 (5) 


TABLE IV. 


tt'^ contribution to muon ^ — 2 in the VMD b model. 




r(7r° ^ 77) [eV] 


A^o [GeV] 


a^(7rO)/(a/7r)^ 


7.22 


0.77 


-0.011 86 (2) 


7.22 


0.80 


-0.043 78 (2) 


7.22 


0.82 


-0.045 67 (2) 


7.78 


0.75 


-0.043 83 (2) 


7.78 


0.78 


-0.045 85 (2) 


7.78 


0.80 


-0.047 83 (2) 


8.34 


0.73 


-0.045 75 (2) 


8.34 


0.76 


-0.047 86 (2) 


8.34 


0.79 


-0.049 93 (2) 



29 



TABLE V. rj contribution to muon g — 2 in the VMD b model. 



T{ri^ 77) [keV] 


[GeV] 




0.42 


0.78 


-0.009 645 (4) 


0.42 


0.81 


-0.010 291 (4) 


0.42 


0.84 


-0.010 940 (4) 


0.46 


0.76 


-0.010 066 (4) 


0.46 


0.79 


-0.010 742 (4) 


0.46 


0.81 


-0.011 419 (5) 


0.50 


0.74 


-0.010 467 (4) 


0.50 


0.77 


-0.011 171 (4) 


0.50 


0.79 


-0.011 876 (5) 


TABLE VL 


rf contribTition to nmon g - 


- 2 in the VMD« model. 


r(??' ^ 77) [keV] 


A^, [GeV] 




4.07 


0.82 


-0.008 655 (4) 


4.07 


0.85 


-0.009 295 (4) 


4.07 


0.88 


-0.009 942 (5) 


4.26 


0.81 


-0.008 823 (4) 


4.26 


0.84 


-0.009 475 (4) 


4.26 


0.87 


-0.010 135 (5) 


4.45 


0.80 


-0.008 986 (4) 


4.45 


0.83 


-0.009 651 (5) 


4.45 


0.86 


-0.010 323 (5) 
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TABLE VII. Values of and /? in the extended model (I) for 7r° for the various values of Mx 
and r(7r° 77). The error accompanying each parameter is obtained with the other parameter 
fixed to the best fitted value and at la confidence level. 



Mx [GeV] 


r(7r° ^ 77) [eV] 




/3(7rO) 


Mp 


7.22 


0.12t°;}6 


-0 0003+° °^^^ 

U.UUUO_o Q218 


Mp 


7.78 


o.o2l°:}g 


U.UUOO_Q Q2Q4 


Mp 


8.34 


-0.061}:°^ 


U.UUDD_o.oi94 


2.0 


7.22 


'-'■-'^^-0.12 


-0.0005lHiiS 


2.0 


7.78 


"-'■"-'^-0.10 


0.0024t°:°ii 


2.0 


8.34 


-0 06+°-^° 


OO'^'S+o oise 

U.UUOO_o 0137 


4.0 


7.22 


„ ^ „_i_0 16 

o-i2io:i2 


n nnn/l+0 0121 

-0.0004lo-u^^^ 


4.0 


7.78 


02+°-^° 


o.oo2ol°:°ili 


4.0 


8.34 


-0 0fi+°-^2 


0.0040lHill 


6.0 


7.22 


"-'■-'^^-0.12 


-0.0003lHil^ 


6.0 


7.78 




o.ooi8t°:°igl 


6.0 


8.34 


-0 06+°-^° 


0.0036inig^ 
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TABLE VIII. Values of and (3 in the extended model (I) for rj. The meaning of numbers 
are the same as in Table VII. 



Mx [GeV] 


r(r/ — > 77) [keV] 






Mp 


0.42 


-0 12+°-^^ 


n5q+°-°22 

u.uoy_Q 


Mp 


0.46 




U.UOU_Q Q2Q 


Mp 


0.50 


-0 25+°-°^ 


U.UOi_o.019 


2.0 


0.42 


_0 iq+0.09 


U.UOJ_o.015 


2.0 


0.46 


-0 25+°-°'^ 


0.05410:°!^ 


2.0 


0.50 


-o.29:Hi 


0.054l°:°Jl 


4.0 


0.42 




n n4n+°-°^3 

U.U4U_o.013 


4.0 


0.46 


fi oo+U.Uo 

'^•^•-'-0.06 


U.U4U_o.oi2 


4.0 


0.50 


n 97+0.06 
-0.05 




6.0 


0.42 


-0 15+°-^° 
'^•-'^'^-0.08 


U.UOO_o.012 


6.0 


0.46 


-0 21+°-°^ 


U.U,50_Q Q^]^ 


6.0 


0.50 


-n 26+°°^ 

^•^"-0.06 


U.UOD_o oil 
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TABLE IX. Values of and (3 in the extended model (I) for Vj' . The meaning of numbers are 
the same as in Table VII. 



Mx [GeV] 


r(V ^ 77) [keV] 






Mp 


4.07 


•J-jU-o.is 


n ni5+°-°^'^ 

U.UiO_Q Q^7 


Mp 


4.26 


n 00+0.22 
'-'•^^-0.16 


U.UiO_o.oi6 


Mp 


4.45 


n 1 5+°-i9 


Q Qi 7+0.016 
'-0.016 


2.0 


4.07 


n 25+0-23 


n 01 f^+o-013 


2.0 


4.26 




U.UiO_o.oi3 


2.0 


4.45 




U.UiD_Q Q^2 


4.0 


4.07 


n 25+°-23 


012+°-°^° 

U.UiZ_Q Q^Q 


4.0 


4.26 


n 1 7+0.20 


n ni 2+0-010 


4.0 


4.45 


11+0-17 


Q QlO + O.OlO 


6.0 


4.07 


-0.17 


nin+°-°°^ 

U.UiU_Q QQ9 


6.0 


4.26 


1Q+0-20 


f.-|f.+0.009 
U.UiU_o.009 


6.0 


4.45 


n is+°-^s 

'-'•l'J-0.13 


Q +0.009 
U.Uii_o.009 
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TABLE X. 7r° pole contribution to muon ^ — 2 in the extended model (I) . The four numerals 
corresponding to one set of k^) are the values of a^/{lO~'^{a/'JT)^) for the central value of 
T{7r^ — > 77), and, from the above, for Mx = Mp, 2.0, 4.0 and 6.0 GeV respectively. The errors 



accompanying them are inferred from the uncertainty in r(7r0 — > 77). 





/3(7r°)min 


/3(7r°)center 








A qq-, +0.203 








-4 402+°-2i4 

4.4UZ_Q 220 

-4 407+0-220 

-4 41 1+0-225 
4.411_Q 225 








4.00y_Q 208 


-4 621+°-227 

4.DZi_Q 231 


I^(j>{j^ ) center 


-4 442+0-186 

-4 417+0-191 
^•^-■-'-0.199 


-4 fi70+°-2i2 

-4 t^7fi+0-216 


-4.70ll°J|° 
-4.736lg:ii 




-4 400+°-^^^ 

4.4UU_Q 202 


-4.579l0;i3 


-4.761^°:i? 






4 70-1+0.201 

'I. (0±_Q 208 




'^0(7I"'^)max 




-4 742+0-211 

^■'^^-0.216 

-4.748l°:i^ 
-4.7511°:!^ 
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TABLE XI. r] contribution to muon g—2 in the extended model (I). The numbers corresponding 



to one set of (/?, k^) have the meaning similar to those of Table 



/3(?/)mm 


/3(??)center 


/5(??)max 




-n QQ^+0 063 
U.jyD_o.o62 




^^(^)min 


i.±iy_Q Qg3 
_1 1^^4+0.085 

1.104l_Q Qgg 

-1 1 74+0-087 
-'-•-'^'*-0.087 




-0.990l°:°i 


1 n-^il +0.063 
i.UOi_o.o62 


-1 074+0-071 

-•-■U' "±-0.070 


^(f>\V jcenter 

-1 090+°-°'^^ 


-1 146+0083 

1 100+0.086 
J--J-O^_0.085 


1 99Q+0.094 
i.ZZy_Q 

1 978+0.098 
-"-•^'"-0.098 


i.iU4_g Q7g 


1 009+0.087 

i.ZUZ_Q Qgg 


1 *^ns+o-ioi 

-i.rfU»_o.lOO 




-1 069+0064 




'-5'(/)(^)max 


1 170+0.083 
-'■•^'•-'-0.084 

-, 9ir,+0-086 
l.ziu_Q og7 

-1 9qn+0.088 
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TABLE XII. r/' contribution to muon g' — 2 in the extended model (I). The numbers corre- 



sponding to one set of Kfp) have the meaning similar to those of Table 



/3(??')mm 


PiV )center 


/3(V)max 




-n 842+°-0i^ 




K(k(71 )min 


-0.859lo:o2i 
-0.8721°:°^} 
-0.87910:01 






-0.89710:018 


-0 906+0-021 

U.yUD_o 021 


^<f>\V /center 


-0.910lO:0|i 

-0 Q2S+0-022 


-0.94210:01^ 
-0.96510:026 


-0.88Atro\l 


-0 qsO+0-022 

U.yOU_Q Q22 


-0.97710:02V 




u.yoo_o 019 




l^ipiV )max 


-0.96210:021 

U.y (0_Q 022 

-0 982+0-023 
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TABLE XIII. f3 and in the extended model (II) for tt" for the various values of r(7r'' — > 77). 



The error accompanying each parameter 


is obtained with the other parameter 


fixed to the best 


fitted value and at la confidence level. 






r(7r0 ^ 77) [eV] 




/3(vrO) 


7.22 


^•J^"J-0.12 


-0 004+0-042 

U.UU4_Q Q43 


7.78 


'-'•'-''-0.10 


n nil +0.043 
U.Uii^Q Q43 



8.34 -0.091°:°^ 0.0281^41 



TABLE XIV. /5 and in the extended model (II) for r]. The meaning of numbers is the same 



as in Table XIII 



r(r/ ^ 77) [keV] 






0.42 


n 90+0.O6 


SO+0 08 
'-'•'J '-'-0.07 


0.46 


-0 SO+0-05 

'-'••J '-'-0.05 


n 97+0.08 
'-'■-^' -0.07 



0.50 



-0.32 



+0.05 
-0.04 



0.24 



-0.07 
-0.06 



TABLE XV. 13 and 
as in Table XIII. 


in the extended model (II) for ij' . 


The meaning of numbers is the same 


r{r]' 77) [keV] 




/3(r/') 


4.07 


-0 ■^1+0-05 
'-'•"5-1-0.05 


20+0-03 

'-'•^'-'-0.03 


4.26 


_n Qn+o-05 

'-'••3U-0.05 


n 1 Q+0-03 
^■1-^-0.03 


4.45 


n 90+0.05 


1 7+0-03 

'-'•1 ' -0.03 
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TABLE XVI. TT^ pole contribution to a^/(10 ^(a/Tr)^) in the extended model (II). The errors 
accompanied with them are inferred from the uncertainty in r(7r'^ — ^ 77). 





/3(7r'^)min 


^(71"°) center 


/3(7r°)max 


H(l>\'^ Jmin 




-4.386l°:Ji 




"'9V" ^center 


-0.189 


-4.559l°:i^^ 


-4 fifi2+°-i96 


Juiax 




-4.736t°:Sl 




TABLE XVII. 


77 pole contribution to a^/(10 ^(ck/tt)^) in the extended model (II). 




/3('?)min 


/?('?)center 


/?(??)max 






U.05Z_Q Q4g 




'*'0V'// center 


U-"' ' -0.052 


-0.6381°:°^^ 


-0.6001°:°^? 






-0.697lHii 




TABLE XVIII. 


f]' pole contribution to 


a^/(10 ^(a/7r)^) in the extended model (II). 




/3(?7')min 


/3(?7')center 


/?(V)max 


^■(.>(jl )miii 




n 09 1+0.042 




^^(pijl )center 


-0.844l°:°?« 


n 879+0.038 


n SQ'^+0 045 
u.o»o_o 038 


'^0(V)max 




n Q1 c;+0.036 
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FIGURES 

FIG. 1. Leading diagrams in the chiral perturbation and 1/Nc expansion which contribute 
to the hadronic hght-by-hght scattering effect on the muon anomaly: (a) pion-loop diagram, (b) 
pseudoscalar pole diagram, (c) quark-loop diagram. Solid and wavy lines represent muon and 
photon respectively. Dotted line in (a) corresponds to the charged pseudoscalar meson while 
dotted line connecting the two blobs in (b) corresponds to the neutral one. The closed solid line 
in (c) represents the quark loop. These diagrams are typical ones in respective classes. 

FIG. 2. Pseudoscalar pole diagrams. Dotted line between the two blobs represents the prop- 
agation of neutral pseudoscalar meson. 

FIG. 3. Comparison of various theoretical form factors of r(7r'' — > 77*) with the experimental 
data. Solid line corresponds to the VMD model (61) while dashed line (62) and dot-dashed lines (63) 
correspond to the QTL model and the QTL model supplemented by the VMD model, respectively. 

FIG. 4. Comparison of the various theoretical asymptotic behavior of form factors of 
r(7r'' — > 77*) with the experimental data. Three straight lines, from the above, are the Brod- 
sky-Lepage (BL), OPE and the Gerard-Lahna (GL) predictions, respectively. The curve represent- 
ing the Brodsky-Lepage interpolation formula ( |2.13 ) is also shown for comparison. 
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